Degeneracy of the ground-state of antiferromagnetic spin-1/2 Hamiltonians. 
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In the first part of this paper, the extension of the bieb-Schultz-Mattis theorem to dimensions 
larger than one is discussed. A counter example to the original formulation of bieb-Schultz-Mattis 
and Affleck is exhibited and a more precise statement is formulated. The degeneracy of the ground- 
state in symmetry breaking phases with long-range order is analyzed. The second and third parts 
of the paper concern resonating valence-bond (RVB) spin liquids. In these phases the relationship 
between various authors approaches: baughlin-Oshikawa, Sutherland, Rokhsar and Kivelson, Read 
and Chakraborty and the bieb-Schultz-Mattis- Affleck proposal is studied. The deep physical relation 
between the degeneracy property and the absence of stiffness is explained and illustrated numerically. 
A new conjecture is formed concerning the absolute absence of sensitivity of the spin liquid ground- 
states to any twist of the boundary conditions (thermodynamic limit). In the third part of the paper 
the relations between the quantum numbers of the degenerate multiplets of the spin liquid phases 
are obtained exactly. Their relationship with a topological property of the wave functions of the low 
lying levels of this spin liquid phase is emphasized. In spite of the degeneracy of the ground-state, 
we explain why these phases cannot exhibit spontaneous symmetry breaking. 
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I. INTRODUCTION 



The question of the degeneracy of the ground-state of 
spin-i Hamiltonians has a long history. One part of this 
history goes back to the beginning of the sixties when 
Lieb, Schultz and Mattis proved that a spin-^ antiferro- 
magnetic periodic chain ofJength L has a low energy 
excitation of order 0(1/ Lp. This theorem (called in 
the following LSMA) was then extended by Lieb and 
Affleck to half-integer spin but fails for integer onesH. 
It states that in one dimension (ID), SU(2) invariant 
Hamiltonians with half-integer spins in the unit cell, ei- 
ther have gapless excitations or degenerate ground-states 
in the thermodynamic limit. The next important contri- 
butions to this long debated questioftrCaHie at the end 
of the eighties when various authorscloBQ showed that 
the ground-state and first excitations of short-range Res- 
onating Valence-Bond (RVB) spin liquids with one spin-i 
per unit cell were 4-fold degenerate on a two dimensional 
(2D) torus (periodic boundary conditions). The connec- 
tion between the two statements in ID was promptly es- 
tablished! In fact Lieb Schultz and Mattis had suggested 
in their original paper that their ID theorem could prob- 
ably be extended to higher space dimensions: to support 
this conjecture, they developed an heuristic argument, 
that has been refined later on by Afflecktl. 

In this paper we examine the various arguments that 
have been formulated for a 2D spin systems at the light 
of a large number of numerical results. 

In section II, we study the situations with long-range 
order (either magnetic Neel long-ranged order or Valence 



Bond Crystals) and we show that these systems indeed 
have a degenerate ground-state in the thermodynamic 
limit. We study two examples in which one aspect of the 
LSMA conjecture appears to be incorrect in 2D. 



In section |HI|, we go.through the analysis of a recent ar- 
gument by OshikawaH We present a numerical analysis 
illustrating the relative part of exact results and conjec- 
ture in Oshikawa's argument and we show numerically 
that for the MSE spin liquid Oshikawa's assumption is 
verified and the ground-state is 4-fold degenerate. 



In section IV we consider the general case of short- 
range RVB spin liquids (RVBSL). In such systems the 
degeneracy is intimately related to the existence of topo- 
logical sectors in the Hilbert space of short-range dimcr 
coverings. In such a framework we determine the re- 
lations between the quantum numbers of the degenerate 
multiplets both for even x odd samples and for even x even 
ones. These results are illustrated on the quantum hard- 
core dimer (QHCD) model on the triangular lattice. We 
finally explain why in spin liquids spontaneous symmetry 
breaking is impossible in spite of the degeneracy of the 
ground-state. 



In the last section (section |v|) , we summarize the most 
salient results of the present work and conclude. Some 
technical points are presented in two appendices. 



2 



II. THE LSMA CONJECTURE FOR 2D 
SYSTEMS. 

A. The LSMA theorem for chains and stripes of 
vanishing aspect ratio. 

The LSMA theorem applies to spectra of quantum 
antiferromagnetic chains with periodic boundary condi- 
tions, and half-integer spin in the unit cell. The proof 
of the theorem is based on the construction of a low ly- 
ing excitation, which may be pictured as a slow twist of 
the exact ground-state. Precisely the excitation is deter- 
mined by the action on the exact ground-state \4>o) of the 
unitary operator U defined as: 



U((j>) = exp(i 



L 



E E< 

j=0 71=0 



qz \ 



(1) 



where j represents the z component of the spin opera- 
tor at site (n, j), the chain being of length L x and width 
L y . The low energy variational excitation reads: 



gLSMA 
7 2tt 



) = U(2tt) |Vo) • 



(2) 



We will design as a column of the sample all the spins 
swept by the unit cell subjected to the L y translations 
along y. Let us call T x the operator for one-step trans- 
lation in the x direction: the Hamiltonian is supposed 
to be translationally invariant [H, T x ] = 0. It is easy to 
show that U (2tt) anti-commutes with T x as flPjOn as the 
number of spins in a column is an odd integertlll. In these 
conditions, if the ground-state wave- vector is ko, the vari- 
ational excitation |#2f MA ) nas a wave- vector ko + (tt, 0): 
it is thus orthogonal to the ground-state. 
The energy of this variational state is: 
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where Eq = (ipo\ H.q \ipa) is the ground-state energy and 
a a finite quantity of order 0(1) measuring the average 
value in the ground-state of spin-spin correlation func- 
tions. The LSMA theorem concerns short range Hamil- 
tonians: for the simple nearest neighbor Heisenberg case 
a is equal to —2/3 times the energy per bond in the 
ground-state. 

Eq. (||) expresses the LSMA theorems: if the system 
under consideration is a chain, \&2-k MA ) ^ s orthogonal to 
the ground-state and collapses onto it as 0(1/ L x ) when 
the chain length L x goes to infinity. The property re- 
mains true for a stripe with an odd number of rows and 
an odd number of spins | in the unit cell as soon as 
the aspect ratio L y /L x goes to zero when the size of the 
sample goes to infinity: the energy of the excited state 
is then 0(L y /L x ) above the ground-state. On the other 
hand, if the thermodynamic limit is taken with a non 
vanishing aspect ratio, the gap between the ground-state 



and the variational excitation does not close in 2D and 
higher dimensions. n 

Faced to this alternative, Lieb, Schultz and MattisEJ 
expressed the conjecture that the degeneracy observed 
in the odd stripes with vanishing aspect ratio could be 
a true physical property of the 2D samples with aspect 
ratio equal to 1. They wrote: "Because the excitation 
energy of exact low-lying states should not depend on 
the shape of the entire lattice, there should be no energy 
gap for a lattice of N x N sites either". Elaborating 
on this consideration, Afflecktj expressed essentially the 
same conjecture. 

In the following subsections, we want to precise these 
statements at the light of a numerical analysis of the 
exact spectra of various generic 2D quantum antiferro- 
magnets. To discuss the validity of the above-mentioned 
conjecture for dimensions larger or equal to 2, we extract 
from it different proposals with increasing specificity: 

• proposal A: The ground-state of an antiferromag- 
netic system with an odd number of spins i in the 
unit cell is degenerate in the thermodynamic limit. 

• proposal B: If ko is the wave- vector of the ground- 
state there are at least 2 additional eigenstates with 
wave vectors k + (tt, 0) and k + (0, tt) that collapse 
onto the ground-state in the thermodynamic limit. 

• proposal B : In the thermodynamic limit the unit 
cell of the ground-state is enlarged by a factor 2 d , 
where d is the dimension of the lattice. This state- 
ment is a bit more precise that statement B: in 2D 
it implies the degeneracy of the ground-state with 
eigenstates of wave vectors: k + (tt, 0), (0,tt) and 
k + (tt, tt). 

• proposal C: A good variational estimate of the ex- 
cited states of proposal B can be obtained by a 
twis t of t he ground-state (to be specified in subsec- 
tion |n^). 



B. 



Proposal A and spectra of different generic 2D 
quantum antiferromagnets. 



At T = the ground-state of an antiferromagnetic sys- 
tem can be Neel ordered (the spin-spin correlations ex- 
hibit long-range order), a Valence-Bond-Crystal (VBC) 
(the dimer-dimer correlations show long-ij&tige order) or 
an RVBSL state with no long-range ordenlj. 

In any of these cases and for various reasons the 
ground-state is degenerate in the thermodynamic limit. 

In the first two-cases the symmetry breaking of the 
macroscopic ground-state implies a degeneracy of the 
ground-state in the thermodynamic limit. 

i) In the case of Neel long-range order, there are both 
degenerate states (which form the true thermodynamic 
ground-state) and gapless excitations. The gapless ex- 
citations, the antiferromagnetic magnons, are the Gold- 
stone mode associated to the broken continuous SU(2) 
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symmetry. Whereas it is well known that the soft- 
est magnons scale as 0(l/L), it is sometimes taken for 
granted that these collective "first excitations" are the 
first excited levels of the multi-particle spectra. This is 
indeed false: the "T=0 Neel ground-state" (or "vacuum 
of excitations" ) is itself a linear superposition of ~ N a 
eigenstates of the SU (2) invariant Hamiltonian which col- 
lapse to the ground-state as 0(1/N) (a is the number of 
sub-lattices of the Neel state, L is the linear size of the 
sample and N the number of sites of the sample): this 
is the tower oL states, invoked in P. W. Anderson 1952 
famous paperOli3£3liil. 

ii ) In the case of a discrete broken symmetry, as for ex- 
ample the space symmetry breaking associate d .t o . toftSr 
range dimer or plaquette order (VBC)^00li3O'L£lEil, 
or in the case of T-symmetry breaking associated to long- 
range order in chirality, the vacuum of excitations is also 
degenerate in the thermodynamic limit. The spectrum 
of a finite-size sample shows a quasi-degeneracy of the 
ground-state. The number of quasi degenerate states is 
finite and a function of the symmetry of the order param- 
eter, it is independent of the sample size. The collapse 
of these quasi-degenerate levels on the ground-state is 
supposed to be exponential with the size of the lattice. 

in) In the last case of an RVB state with no long-range 
order we have found two generic situations: 

a) A spin liquid with a gap to all excitations with 
up to now two examples: a phase of the multiple-spin 
exchange Hamiltonian on the triangular lattice£2l, and a 
phase pf the frustrated J\ — J2 model on the honeycomb 
latticeE3 (we shall call this phase a type I spin liquid). 
Finite size studies of these cases point to a degenerate 
ground-state on the triangular lattice and a unique one 
on the honeycomb lattice. As explained in section IV, 



this difference comes from the fact that the triangular 
lattice has one spin ~ per unit cell and the honeycomb 
lattice has 2. 

(3) A spin liquid with a spin gapihled with a continuum 
of singlets on the kagome lattices, and on the triangular 
latticeE.il. In the last two cases the continuum of singlets 
appears to be adjacent to the ground-state (type II spin 
liquid). 

In any of these cases Proposal A appears to be true, 
and the restriction to problems with half-odd integer 
spins in the unit cell necessary (as can be seen on the 
example of the J\ — J2 model on the honeycomb lat- 
ticcM). 



C. Proposal B and B : a thorough study of one 
counter-example 

Proposal B is valid when the macroscopic ground- 
state unit cell is multiplied by a factor 2 d . This is for 
example the case of the 4-spin S — plaquette orderj-of 
the J\ — J2 model on the square lattice for J2/J1 ~ 0.5E1 
Proposal B' is not valid in the ground-state of the J%l J\ 
model for J%j J\ » 0.5. In that case the ground-state is 




FIG. 1: Brillouin zone of the triangular lattice: wave- vectors 
Ki = (-7T, 0) frequently invoked in the text are shown in the 
Brillouin zone. When the Cev symmetry of the sample is 
preserved, eigenstates with wave- vectors Ki, K 2 , K 3 are de- 
generate. The two kinds of samples used in subsection 2.3 are 
designed as A, and B. 



a two-sublattice columnar antifcrromagnet which sponta- 
neously breaks the translation and 7r/2 rotation symme- 
tries. It is easy to check that proposal B' is not satisfied 
in that case since the spatial symmetry breaking in that 
phase is incompatible with momentum k = (7r,7r). The 
four-fold ground-state degeneracy is obtained by a com- 
bination of two states with zero momentum, one with 
k = (0,7r) and one with (tt,0). 

We will now show that even the weaker proposal B can 
be false in true 2D spin systems and that, contrary to 
LSMA assumption, a 2D sample with vanishing aspect 
ratio can have a mathematical spectrum different from 
that of the true 2D samples. To illustrate this assumption 
we have done a thorough numerical study of the exact 
spectra of the Heisenberg model on the triangular lattice 
for various shapes and sizes of the samples. 

The Heisenberg Hamiltonian with nearest neighbor in- 
teractions reads: 



Ho — 2 S,; • Sj. 



(4) 



<i,3> 



The Bravais lattice is defined by the two unit vectors 
u and v . The reciprocal lattice is hexagonal: the lattice 
and Brillouin zone are displayed in Fig. ^. The sample is 
defined by two vectors: 



Ti( 2 ) = Zi(2)U + n»i(2)V. 
Periodic boundary conditions read: 



S(R 4 



S(Ri) 



(5) 



(6) 



For a given sample the wave- vectors compatible with the 
boundary conditions depend on Ti,T2. 
pPrevious studies on samples with aspect ratio equal to 
fli-3 have established that this model has a three sublat- 
tice Neel order. The Anderson tower of states collapsing 
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to the absolute ground-state as 1/N encompass eigen- 
states with wave- vectors: k = and 
k = ±(27r/3, 27t/3). The low-lying magnons have wave 
vectors in the neighborhood of the points T and Wi of 
the Brillouin zone of the system (see Fig. |l|). Due to 
the overall consistency of the picture, we suspected that 
eigenstates with wave- vectors (0, it) or (71", 0) (points Ki of 
Fig. |l|) appear only in components of energetic magnons. 

In this work we report a thorough numerical study of 
the behavior of the first k = k + (ir, 0) exact eigenstate 
of this model as a function of the aspect ratio and size of 
the samples (ko is always the wave vector of the absolute 
ground-state for the sample under consideration). 

We have studied samples with various values of 
(L x ,L y ) (even, odd), (even, even) accommodating eigen- 
states with wave- vector (it, 0). For a given aspect ratio 
two kinds of samples have been studied: 

• The simplest one has its sides parallel to the prim- 
itive vectors of the lattice: Ti = L x u, T2 = L y v. 
If L x or L y = 1,2 mod 3 , PBC on such samples 
frustrate the natural Neel order of the lattice: the 
ground-state is usually not ak = state and the 
spectrum is not typical of Neel order. We neverthe- 
less study these samples which are perfectly reason- 
able samples with respect to the LSMA conjecture 
(A sample of Fig. [I]). 

• If one dimension only is not multiple of 3, one may 
use a different shape, with the same aspect ratio, 
where PBC do not frustrate Neel order (see B sam- 
ple of Fig. [l| and double starred samples in Table |). 
On these non frustrating samples (NFS) the abso- 
lute ground-state is a k = state and the spectrum 
has the typical structure of a Neel ordered system. 
On such samples one can also build the variational 
state described by Eq. (g). It has a wave- vector 
(71", 0), and its gap to the absolute ground-state is 
iaLy/Lx (to be compared with Eq. (||) for an A 
sample). 

This study leads to the following observations (see Ta- 
ble|: 

• For a given size the spin gap A s has a weak de- 
pendence with the sample shape (it is especially 
weak on NFS). It decreases with the sample size as 
expected for a Neel ordered system (see Fig. |J): 

XN py/N N z 

where x is the homogeneous spin susceptibility, c 
the spin- wave velocity, p the spin stiffness and /? is a 
number of order one. The first triplet excited states 
of NFS have wave-vectors k = ±(27r/3, — 27r/3) and 
k = 0. 

• The gap Ao to the first excited state with total 
spin 5 = and wave-vector k + (n,0) behaves 
differently. (This state is the first exact state with 
the same wave- vector as |#2,f M ' 4 ))- 
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Lx 


Ly 


A s 


Ao 


Ai 


5 




12* 


Vl2 712 


1.726 


2.471 


2.836 


2.47 


no 


18* 


6 


3 


1.131 


1.483 


2.882 


2.96 


no 


24 


8 


3 


0.942 


0.474 


1.486 


1.26 


yes 


24 


6 


4 


1.193 


1.424 


1.697 


2.13 


no 


2^.** 


8 


3 


0.926 


1.034 


2.073 


2.75 


no 


24** 


6 


4 


0.922 


1.061 


2.175 


1.59 


no 


30 


10 


3 


0.799 


0.344 


1.183 


1.15 


yes 


30 


6 


5 


0.825 


0.762 


1.652 


0.91 


yes 


30** 


10 


3 


0.782 


0.732 


1.618 


2.43 


yes 


30** 


6 


5 


0.776 


1.067 


2.003 


1.28 


no 


36* 


12 


3 


0.692 


0.548 


1.324 


2.19 


yes 


36* 


6 


6 


0.688 


1.359 


1.785 


1.36 


no 



TABLE I: Gaps of various excited states (from the absolute 
ground-state): A s , spin gap; Ao, gap between the absolute 
ground-state (with wave vector ko) and the first excited state 
with wave vector ko + (n, 0) in the S — sector; Ai, gap 
between the absolute ground-state and the first state with 
wave vector ko + (n, 0) in the S — 1 sector; S, is the ratio 
of Ao to the aspect ratio of the sample L y /L x . In the last 
column we indicate if [ko + (7T, 0), S = 0) is the first excited 
state of the spectrum. The lines with one or two * indicate 
A or B samples with PBC compatible with 3-sublattice Neel 
long-range order (See Fig. |l]for the double-starred samples). 



1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — r 




FIG. 2: Finite-size scaling of various gaps to excited states. 
Black circles: A s gap to the first S = 1 state (it has 
wave- vector Wi(2) in the Brillouin zone of Fig. fy). Open 
squares: Ao gap to the first \(ir, 0),S = 0} state (points Ki 
in the Brillouin zone). Black diamonds: Ai gap to the first 
10, 0),S = 1} state. 
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• Ao is small for small aspect ratio and for a given 
size is very sensitive to the aspect ratio and to the 
boundary conditions. For a small aspect ratio this 
state is indeed the first excited state of the exact 
spectrum (last column of Table |) . 

• For a given size and aspect ratio, Ao increases when 
going from frustrating periodic boundary condi- 
tions, to boundary conditions compatible with the 
three sublattice Neel order (except for the 6x4 
sample, where A is always large). 

• For a given size, Ao increases with the sample as- 
pect ratio: whatever the size of the sample this gap 
is always of the same order of magnitude as the as- 
pect ratio (6th column of Table |) 

• Let us consider now the two samples with 6x5 
and 6x6 sites which do not frustrate Neel order 
(double starred and starred lines of Table |). When 
going from the N=30 spectrum to iV = 36 all gaps 
decrease, which seems natural, except the Ao gap 
which increases noticeably. This is a second indi- 
cation of the strong dependence of the position of 
this specific level on the sample's aspect ratio. 

• We finally looked at the size and shape effects on 
the gap Ai to the first state with wave-vector (w, 0) 
and total spin 5=1: | (7r, 0), 5=1) (Let us recall 
that |#27r MA ) has not a fixed total spin). For any 
size and shape Ai > Ao. Ai data corresponding to 
NFS with the largest aspect ratios obey the same fi- 
nite size scaling law as the A s spin gap (FigJ|) : this 
is exactly what is expected for a (tt, 0) magnon exci- 
tation (the magnon excitation could be written as a 
twist of a symmetry-breaking ground-state, which 
itself is a coherent superposition of the low lying 
levels of the Anderson tower of states). Extrapo- 
lation to the thermodynamic limit gives a energy 
of 1.17 (the linear spin wave theory gives the value 
2, but all the ED data point to a softening olthe 
spin- wave spectrum by quantum fluctuationso) . 

• It is difficult to assign a value to the thermody- 
namic limit of the gap A to the first | (tt, 0), 5 = 0) 
state for NFS with the largest aspect ratio (see 
Fig.^j) (two samples only, 12 and 36 have aspect 
ratio equal to 1): however examination of all the 
data pleads in favor of a non zero gap. If the semi- 
classical picture is valid, the state \(ir, 0),5 = 0) 
belongs to the tower of states of the one-magnon 
excitation of wave vector (tt, 0), and in the thermo- 
dynamic limit 

lirriN^oo (A ) = Hmjv^oo (Ai) . (8) 

The present data are too scarce to sustain this last 
assertion but they do not go against it. 

This study suggests that: 



i) The thermodynamic limit: 

L x , L y — > oo with L y /L x — + is special, it could depend 
on the number of rows in the sample (see sample 6 x 
4) and is different from the limit: L Xl L y — > oo with 
L y /L x = 0{l). 

ii) The first state \(n,Q),S = 0) which collapses to the 
ground-state in the quasi one-dimensional system is not a 
"natural" excitation of the 2D system: the natural exci- 
tations are the magnons and the magnon with the wave- 
vector (7r, 0) is not a soft mode as soon as the period of 
the symmetry-breaking ground-state is not even (as in 
the case above, where the period is 3). 

Hi) The nearest neighbor Heisenberg model on the tri- 
angular lattice appears to be a serious counter-example 
to proposal B and B for a true 2D system. We suspect 
that the same model on the square lattice would equally 
be a counter-example to these proposals. We have not 
done a full size and shape analysis of this case as we 
have done for the triangular lattice. But we have ver- 
ified that for samples with non zero aspect ratios the 
Anderson tower of states encompasses eigenstates with 
wave-vectors k = (0, 0) and k = (tt, 7r) as expected, and 
the k = fat. 0) and k = (0, tt) are components of energetic 
magnonsEJ. 



D. Proposal C 

Eq. (j^) implies that in a true 2D system it is in gen- 
eral impossible to build a low-lying excitation using the 
LSMA twist operator U(2n). In RefB, Affleck discussed 
the general possibility to build a one-dimensional excita- 
tion in a two-dimensional medium. We do not know how 
to do this in a general case, and as explained in the pre- 
vious subsection, we suspect that this quest is in general 
doomed to failure. 

However, in the specific case of spin liquids, it could 
be possible to build such an excitation, inasmuch as the 
stiffness of a 2D spin liquid vanishes when the size of the 
sample goes to infinity. 

In an ordered system, the stiffness p, which measures 
the response of the sample L x x L y to a twist 4> of the 
boundary condition in the x direction, can be defined as 

E(<t>) = E(0) + ^p(J^) +0{^). (9) 

with E((f>) the total energy of the sample. Far from the 
quantum critical point the stiffness is usually of the order 
of magnitude of the coupling constant (~ 1/9 for the 
Heisenberg model on the triangular lattice). 

Let us now consider the following one-dimensional per- 
turbation of the ground-state with PBC: a constraint on 
line number one is forced with the twist operator 

2 L. v -i 

u 1 (2^) = exp( ? — Y, (10) 
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FIG. 3: 2-torus with one cut A. 

and the system is supposed to relax around this default 
at fixed momentum (Eq. ( |Io| ) insures that the variational 
state has momentum (7r,0)). 

In an ordered system with stiffness p the variational 
energy of such an excitation could be approximated as: 



(11) 



with (3 a number of order one. In a true 2D ordered sys- 
tem the correction to the ground-state energy can never 
go to zero for a non-zero aspect ratio. On the contrary, 
in a spin liquid sample, with linear dimensions larger 
than the correlation length, we conjecture that the sec- 
ond term of the above estimate goes to zero with the 
stiffness. 

This heuristic argument is in agreement with the be- 
havior of the exact spectra described in the next sections. 
The key ingredient to discriminate between the two cases 
is the presence of a finite (or vanishing) stiffness: a fact 
intimately related to the central hypothesis of Oshikawa 
that we will discuss in the following section. 



III. OSHIKAWA'S ARGUMENT 

In a recent paper! Oshikawa discussed the conse- 
quences of the presence of a gap in a quantum many- 
particle system with conserved particle number on a pe- 
riodic lattice in arbitrary dimensions. Inspired by Laugh- 
lin's topological arguments for the quantum Hall effect, 
Oshikawa argues that a finite excitation gap is possible 
only when the particle number per unit cell of the ground- 
state is an integer. The above-mentioned formulation 
translated for a spin system reads: "A spin system can 
have a gap only if the spin in the ground-state unit cell 
is integer" and Oshikawa statement for a spin-i system 
is "a spin-i system can have a gap only if, the transla- 
tional symmetry is spontaneously broken" : the ground- 
state unit cell is thus adequately enlarged to be consistent 
with the previous statement. 

To support his statement Oshikawa developed an argu- 
ment a la Laughlin involving the adiabatic transform of 
the ground-state under the action of a fictitious magnetic 
flux piercing the ring of the sample submitted to periodic 
boundary conditions (Fig. |3|). In this section we illustrate 
this adiabatic transform on a spin system. This trans- 
form is here realized in the spin language by twisting the 



boundary conditions (subsection [II A). We then show 
the general mathematical properties of the spectrum of 
any spin Hamiltonian under such an adiabatic transform 
(subsection 3.2 and Appendix A). 

Oshikawa's result is based on a crucial hypothesis that 
the gap from the ground-state manifold to excited state 
does not c lose during the adiabatic transform. In subsec- 
tion [II C we illustrate this assumption on exact spectra 



of the multiple-spin exchange (MSE) Hamiltonian, in a 
range of parameter where this model exhibits a RVBSL 
phase in the thermodynamic limit (see Appendix B for 
the precise definition of the model and the range of pa- 
rameters where it exhibits a spin liquid phase) . We post- 
pone the discussion of the second issue of Oshikawa's pa- 
per, i.e. the discussion of spontaneous symmetry break- 
ing to the next section. 



A. Adiabatic transform of a spin-system under the 
twist of the boundary conditions. 

If the spin-i degrees of freedom are represented by 
hard-core bosons (Holstein-Primakoff), coupling these 
bosons to a fictitious flux is completely equivalent, in 
the spin language, to a twist of the boundary conditions. 
Generalized twisted boundary conditions are defined by 
the choice of the twist axis (here the z axis in the original 
Bq spin frame) and a set of 2 angles <f>j as: 

S(Rj + Tj) = e l ^ s2! ( R Os(R 1 ) e -' i ^^( R ^ (12) 

(see Fig. [j] for the definitions of T\ and T2). For the sake 
of simplicity we will illustrate the general properties of 
this transform on the case of a twist (f> in the Ti direc- 
tion for the simplest SU(2) invariant Hamiltonian: the 
nearest neighbor Heisenberg Hamiltonian (extension to a 
more complex and relevant Hamiltonian is just a ques- 
tion of notations). With periodic boundary conditions 
the Hamiltonian in the Bq spin frame reads: 

Ly-l L x -1 

p=0 n=0 

(13) 

A twist </> in the Ti direction implies the calculation of 
the eigenstates of 



Ly-l 

E 

p=0 



Sij.-i.p • So p + S Lx _i iP • S Li> _x,p+i[L v ]) 



+ 



L y -1 L x -2 

E E 

Ly-l 



where 



p=0 



(14) 



(15) 
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The 4> twist term in Eq. (|4|) is equivalent to the phase 
factor induced by a magnetic flux on the kinetic term 
of charged particles in Laughlin's transform. Under an 
adiabatic twist of the boundary conditions the spectrum 
of 7^0 (Eq. (|T4|)) evolves periodically with a period 2ir. 
But the eigenstates evolution could be more complicated. 
There is no guaranty that the many-body ground-state 
of 7i<p = o adiabatically transforms into the ground-state 
of Ti.tj,=2TT- As we will show below, this is generally not 
the case and the true periods of the eigenstates are in 
general multiples of 2ir. 



B. Twisted boundary conditions and translation 
invariance 

To follow adiabatically an eigenstate of a many particle 
spectrum while twisting the boundary conditions could 
seem difficult because of level crossings during the twist. 
In fact there is an unique, unambiguous way to do it as 
these eigenstates belong to one-dimensional representa- 
tions of a translation symmetry group. Diagonalization 
in such subspaces leads to non degenerate eigenstates 
which never cross. The translational invariance of the 
problem is hidden in its representation in the Bo frame 
but can be restored by a gauge transform. 

Let us define B^, the new frame deduced from the orig- 
inal Bo by the spatially dependent twist described by the 
unitary operator: 

U(<f>) = exp(z-^- £ ]T nS z ntP ) (16) 

x p=0 n=0 

In this new frame, the twisted Heisenberg Hamiltonian 
reads: 



(17) 



(in this equation and in the following, we put a tilde 
on the observables measured in the B^ frame. We will 
indicate by a superscript <f> the eigenstates of 7i^). This 
unitary transform is chosen so that the boundary terms 
in Eq. (04) disappear: 



, Ly-l L x -1 

p=0 n=0 



Tick is translation invariant: 



(18) 

Ti 4,, %( y ) = 0, where T x(y) 
are the operators for one-step translations in the x (resp. 
y) directions. Its spectrum is indeed identical to the spec- 
trum of H.<p. After this gauge transform, we can now 
define one-dimensional irreducible representations of the 
translation group labeled by their wave-vectors k in the 
Brj, frame and adiabatically follow a k eigenstate of Ti.^ in 
the successive Bj, frames while increasing <f>. This is what 
is done in Figs. |5] and ^[ where we see that the true period 
of the many-body ground-state |k = (0,0),5f ot = >^ 



10.5 



10 



Heisenberg Triangular J=l N = 24 



9.5 



8.5 




0.2 0.4 0.6 0. 
0/(2tt) 



FIG. 4: Low-lying spectrum of the 8x3 sample of the Neel or- 
dered Heisenberg model on the triangular lattice as a function 
of the twist <f> at the x boundary of the sample. Circles are for 
\Stot = 0,k = >^ states and squares stand for states with 
\Stot = 0, Ki >^ (wave-vector (ir,0)). Triangles represent the 
lowest Stot = 1 state (whatever its wave vector). Continuous 
lines are drawn through the first \S'tot = 0, k = >^ state, 
the first \Stot = 0, Ki >^ and the lowest Stot = 1 state thus 
giving the spin gap. Due to the finite size, the k wave vector 
of this lowest Sf ot = 1 state changes with the twist angles 
(explaining the discontinuities in the slope of the curve) and 
with the sample size. 



and of the first exact |k = (tt, 0), S* ot — >^ state of the 
8x3 sample (resp. of the 4x5 one) is 4-7T. 

This property is an exact property of the 5 t z ot = 
eigenstates on even x odd samples. A straightforward 
algebra (computation of T x ( y )V {2tt)T^\) , shows that 
U(2tt) maps a state with wave-vector k and total spin 
S^ ot onto a state with the same ^-component of the total 
spin but with a wave vector k' given byc3: 



ki, 



k x 
kit . 



27r^2i + 7r [L y mod 2] 



(19) 



From these relations one immediately deduces that 
a 2tt twist on an even x odd sample maps the absolute 
ground-state |k = (0,0), Sf ot — > of Ho onto the first 
|k = (tt, 0), S% ot = > eigenstate and reversely, whatever 
the nature of the Hamiltonian and the intrinsic proper- 
ties of its ground- state. Fig. ^ is a spectrum of a Neel 
ordered phase and Figs. || and || of spin liquid phases. 

More generally the S% ot = eigenstates have periodic- 
ity 47T on even x odd samples (see Figs. |[ || and ||) and 
2-7T on even x even ones (Figs. [| |?]) . Let us finally remark 
that all eigenstates of a given problem do not have the 
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FIG. 5: Low-lying spectrum of the 4x5 sample of the 
multiple-spin exchange model (Eq. ([is])) as a function of the 
twist 4> at the x boundary of the sample. Same symbols as 
in Fig. We suspect that for this spin liquid, in the ther- 
modynamic limit the wave vector of the first S'ot = 1 state is 
(tt/2, 0). But in the finite small samples that are available we 
never have this wave vector (except in the 4x6 sample): this 
explains the discontinuities in the shape of the lowest S* ot = 1 
state. 



MSE Triangular J a =-2 J 4 =l N=24=4x6 




0/(277) 



FIG. 6: Low-lying spectrum of the 4x6 sample of the 
multiple-spin exchange model (Eq. (E^) as a function of the 
twist (j> at the x boundary of the sample. Same symbols as in 
Fig|. 



FIG. 7: Low-lying spectrum of the 6x4 sample of the 
multiple-spin exchange model (Eq. ([sif)) as a function of the 
twist <j> at the x boundary of the sample. Same symbols as in 
Figg 

same periodicity: the periodicity depends both on the z 
component of the total spin of the eigenstate under con- 
sideration and on its wave-vector (Eqs. (|Tg|) ) . Additional 
information on the quantum numbers associated with the 
point-group symmetries can be found in appendix A. 



C. Robustness of the gap in presence of twisted 
boundary conditions in spin liquids 

To establish the existence of a ground-state degener- 
acy in a gapped system, Oshikawa makes the following 
assumption: the ground-state of Ti. can only be trans- 
formed in itself or into another degenerate ground-state. 
In other words, he assumes that the gap from the ground- 
state manifold to excited states does not close during the 
adiabatic evolution. Thanks to the results of the pre- 
vious sub-section, this hypothesis directly implies that 
in evenxodd samples the \Sf ot = 0,k = (0,0) > abso- 
lute ground-state is degenerate with the lowest \Sj; ot = 
0, k = (tt, 0) > eigenstate. // we take for granted that 
in the thermodynamic limit the spectra of evenx odd and 
evenxeven samplesMre identical, Oshikawa's assumption 
implies Proposal Bo. 

Numerical results for the RVBSL phase of the MSE 
model support these two hypotheses (see Appendix B for 
the precise description of the model and spectral results 
in Figs. § | 0, I and Table ||. The gap A between 
the absolute ground-state \S^ ot = 0, k = (0, 0) > and the 
first \Sf ot = 0, k = (tt, 0) > decreases steadily with the 
size of the sample with no difference between evenxodd 
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MSE Triangular J 3 =-2 J =1 N = 30 



116.5 



117 



-117.5 



118 



I = 1 = 1 B * l-4-i-4-l * 2 





0.4 0.6 
0/(2tt) 



FIG. 8: Low- lying spectrum of the 6x5 sample of the 
multiple-spin exchange model (Eq. ([ii])) as a function of the 
twist 4> at the a; boundary of the sample. Same symbols as 
in Fig. [l|. On this small sample the lowest Stot — 1 state is 
a \S Z = l,k = (27r/3, 0) >. We suspect that in the thermo- 
dynamic limit the lowest S z = 1 eigenstate has a wave-vector 
belonging to the star of (n/2,0) (these wave-vectors are not 
allowed in the present sample). 



N 


20 


24 


28 


30 


36 


A 


4.0027 


2.2944 


1.5683 


0.7026 


0.2527 


A 


0.9967 


0.8195 


0.6958 


0.5615 




A s 


2.9168 


2.1452 


1.2555 


1.1607 


0.9847 



TABLE II: Spin liquid spectra (MSE model defined in Ap- 
pendix B) . Ao measures the gap between the absolute ground- 
state with wave vector (0, 0), and the first excited state with 
wave vector (n, 0) in the S = sector; A measures the in- 
crease in the absolute ground-state total energy in a n twist 
of the boundary conditions; As is the spin gap. 



and even x even samples (see Table || and Fig. ||). 

These results pull into light the central conjecture of 
Oshikawa, i.e. what Laughlin calls the robustness of the 
gap or the uncompressibility of the spin liquid and its 
relationship with our conjecture of the vanishing stiff- 
ness of the 2D spin liquids in the thermodynamic limit. 
Strictly speaking numerical calculations do not allow a 
consistent determination of the "stiffness" of a finite size 
sample in a spin liquid phase: depending on the shape of 
the sample and on the direction of the applied twist the 
response of the system can vary in a range going from 
half to one hundredth of the coupling constant (contrary 
to the case of a Neel ordered system where the stiffness 
is very well defined and totally independent on the di- 



4 
3 

<f 2 
1 





~i 1 r 



1'24 



,28 



36 



>Z2 

^30 



0.02 0.04 

1/N 



0.06 



FIG. 9: Ao gap of the first multiplet in the singlet sector 
versus \fN: the line is an exponential fit oc exp(— v~N/£) 



(£ ~ 0.6) 



rection of the twist). Most of the results obtained on 
the 28, 30, 32-sites samples are 10 to 100 times smaller 
than those of the Neel ordered phase. But the conjecture 
that can be formed from the above mentioned results 
is even stronger than the vanishing of the stiffness and 
amounts to an absolute absence of sensitivity of the ther- 
modynamic ground-state of a spin liquid to any twist of 
the boundary conditions (i.e. all terms of the Taylor ex- 
pansion of Eq. (JsJ) are zero in the thermodynamic limit, 
not only the </> 2 term). 

This conjecture is based on the following observations 
concerning spectra of the MSE model (see Table || and 
Figs 1 ill): 

• in even x odd samples the ground-state energy of 
the system under a twist oscillate between the en- 
ergy of the absolute ground-state at zero twist and 
that of the first excited-state \Sf ot — 0, k = (n, 0) > 
reached for a 27r twistH. The amplitude of the os- 
cillation measured by the energy gap Ao between 
these two states decreases steadily with the system 
size. 

• In even x even samples the amplitude of oscillation 
is reduced: the maximum of energy of the ground- 
state is obtained for a n twist. The amplitude of 
this oscillation A' Q given by the difference in total 
energy between these two states is given in Table ||. 
As the previous gap, this amplitude decreases with 
the system size. 



• The number of samples is too small to ascertain an 
exponential decrease to zero with system length, 
but this conjecture remains plausible (Fig. [)|). 

On physical grounds, the absence of any response to 
a twist is a rather natural property in an RVBSL since 
there such a phase offer no local order parameter this 
perturbation could couple to. In section IV we will come 



back to this question from the point of view of the short- 
range RVB picture of spin liquids. 
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D. Oshikawa's conjecture and Valence Bond 
Crystals 

The numerical work presented in the previous subsec- 
tion concerns an RVBSL. We have not investigated the 
behavior of a VBC with twisted boundary conditions but 
it is clear that the stiffness should be zero. In principle 
non-linear effects could be different in VBC and RVBSL. 



early independent vectors: an exact demonstration exists 
for the coverings of the square latticetil, numerical stud- 
ies have, shown that the property is true for the kagome 
lattice23, and in the triangular lattice for any size large 
enough to insure than any periodic image of a site is at 
a distance larger than 4 units (this work). 

1. Definition of the topological sectors 



IV. THE SHORT-RANGE RVB PICTURE OF 
THE GROUND-STATE MULTIPLICITY OF A 
SPIN LIQUID 

Because of the short-range and antiferromagnetic na- 
ture of magnetic correlations and spin-rotation invari- 
ance, short-range dimer configurations offer a natural de- 
scription of low-energy degrees of freedom in spin liquids. 
Indeed, much insights into their physical properties come 
from the investigation of dimer models. The quantum 
hardj-jfore dimer (QHCD) model of Rokhsar and Kivel- 
sonB'tffl on the square lattice is the simplest dimer model 
but one has to_consider the triangular lattice (Moess- 
ner and Sondho) to find a genuine liquid with a finite 
correlation length. The short-range dimer point of view 
also gave a significant insights on another different kind 
of spin limiid: the Kagome Heisenberg antiferromagnet 
(see Ref. Lj an d references therein). 

First (IV A) we recall that the Hilbert space of short- 
ranged dimer coverings on a torus is made of 4 topolog- 
ical sectors which cannot be coupled by any local oper- 
ator. We show that any wave function described in a 
short-range dimer basis is two-fold degenerate on large 
enough even x odd and can be two- or three- fold degener- 
ate on even x even samples with appropriate point-group 
symmetries. This is illustrated on spectr a of t he QHCD 
model on the triangular lattice. Then (IVB) we argue 
that for a spin liquid where all dimer correlations decay 
exponentially fast with distance this degeneracy is ex- 
tended to 4 for large enough samples. This property is 
checked on the QHCD model and in MSE spectra. 

After some remar ks on the four-fold degener acy in the 
liquid phase ( IV C ) we explain in paragraph [V D why 
RVBSL, in spite of the degeneracy of their ground-states 
cannot manifest spontaneous symmetry breaking. 



A. Decomposition of the space of short-range 
dimer coverings in topological subspaces and exact 
degeneracies 

The space of nearest neighbor dimer coverings of a 
torus is a small fraction of the total 5 = sector 
of the full Hilbert space which increases as a" (with 
a= 1.1753, 2 1 / 3 , 1.3385, 

1.5351 respectively on the hexagonal, kagome, square and 
triangular lattices). These different coverings are not or- 
thogonal, but for large enough sizes they form a set of lin- 



Let us draw a cut A encircling the torus created by 
periodic boundary conditions (see Fig. ||). This hyper- 
surface of dimension d — 1 cuts bonds of the lattice but 
there is no site sitting on it. The position of the cut is 
arbitrary. We may decide in order to follow closely our 
previous discussion to put it between spin L— 1 and spin 
of each row of the lattice. The family of nearest-neighbor 
dimer coverings can be decomposed into two subspaces 
£j± depending on th©,parity Ha of the number of dimers 
crossing the cut Alia. By considering a set of d cuts 
Ai=i,...,d encircling the torus in all possible directions 
one obtains 2 d families of dimer covering. 

Any movement of dimers can be represented as a set of 
closed loops around which dimers are shifted in a cyclic 
way. A local operator will only generate contractible 
loops which will cross each cut a even number of times. 
The number of dimers crossing the cut can therefore only 
be changed by an even integer and the parities Ha; are 
unchanged. 

This property remains true as long as one works in a 
subspace where the dimer lengths are smaller than the 
linear system size, that is when the topological sector 
are well defined (if a dimer length is half the linear of 
the system one cannot decide by which side it goes). 
On the other hand, we checked on the triangular (resp. 
Kagome) lattice that these 4 sectors are the only topolog- 
ical sectors: local 3- (resp 4-) dimer moves can be used 
to transform a configuration of a given sector into any 
other configuration of that sector. 

We show that these subspaces are orthogonal in the 
thermodynamic limitQ. The graph of the scalar product 
< c + |c~ > of two dimer configurations belonging to dif- 
ferent subspaces £ + and £~ has at least one long loop 
encircling the torus in the L x direction. When L x goes to 
infinity this contribution to the scalar product is smaller 
than 2~ La =/ 2 . Consider two normalized vectors \^ + ) and 
l^ - ) belonging to two different sectors: 



/ 



E 



* ± (c ± )|c ± ) 



(20) 



Because of the exponential number of dimer coverings in 
each subspace it is not obvious that and \^~) are 

orthogonal in the thermodynamic limit. It is indeed the 



E 



*+*(c+)*-(c^) < c+\c- > 



(21) 
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< 2- L ' 2 |* + *(c + )*"(c-)| (22) 

c+,c- 

< 2 - i /y^|vi/+( C +)i 2 ^|vi/-( C -)^3) 

< 2" i/2 ||*+|| ■ ||*-|| (24) 

where the norm in the last inequality refers to the diag- 
onal scalar product < Ci||c2 >= S cltC2 . We now wish to 
show that these later bound is finite. In order to see that 
it is the case one can expand the usual scalar product 
< c\\C2 > in powers of x = 2~ 1 / 2 in the spirit of Ref.Q. 
The zeroth order term is nothing but the diagonal scalar 
product: 

<ci||c 2 >=<5 Cl , C2 +0(a; 4 ) (25) 

The convergence of this loop expansion implies that 
||*|| is finite if | (*|*) = 1 and therefore (\fr + |\Ij~) = 
0{2~ L I 2 ). From this we see that not only scalar prod- 
ucts of dimer configurations in different topological sec- 
tors vanish exponentially but this is also true for any pair 
of states. 

In the following, unless explicitly mentioned, we con- 
sider the 2D case for simplicity but most of the topologi- 
cal arguments about dimer covering immediately extend 
to higher dimensions. 

2. Two-fold degeneracy in evenx odd samples 

In the special case of tori with an odd number of rows 
(and an odd number of spin-i per crystallographic unit 
cell), one step translation along the x axis (called T x in 
the following) maps £ + on £~ and reversely. Some point- 
group symmetry can also do this job. A tt rotation about 
a lattice site nearby the cut (called TZ^ in the following) 
has the same effect. If the cut is chosen parallel to a 
symmetry axis of the sample, a reflection with respect to 
this axis (called T. y in the following) will equally map £ + 
on and reversely. 

All these symmetry operations isolate a single column 
C of lattice sites between A and its transform A'. In 
that case columns have an odd number of sites and an 
odd number of dimers must connect some sites inside C 
with sites outside C. Therefore IIa differs from IIa' and 
the two subspaces £ + and £~ arc exchanged. 

For a large enough system these two sectors are 1) 
orthogonal, 2) uncoupled by any local Hamiltonian and 
3) exchanged by symmetry operations (evenx odd). This 
is enough to insures that they have the same spectrum, 
irrespectively of the physics of the model, provided it 
can be described in the short-range dimer space. In fact, 
quantum numbers of these doublets of degenerate states 
are fixed by symmetry. 

We decompose an eigenstate \ipQ > on the two topo- 
logical subspaces defined relatively to the cut A: 

Ho) = + > ( 26 ) 



where belong respectively to the sets \tpo), as 

an eigenstate of the Hamiltonian with periodic boundary 
conditions, belongs to an irreducible representation of the 
translation group. In the following we will also assume 
an Tin and S y invariance of the Hamiltonian and that, 
for simplicity, \ipo > transforms under a one-dimensional 
representation under TZ^ and T, y . 

ft* l^o> - Po Ho) 

W |^o) = ol |^o) (27) 

In the thermodynamic limit |^o") and I^q") are orthogo- 
nal and T Xl TZ^ and S w map £ + on £~ and reversely: 

%\^)=e^\^) 
K„ |Vo=) = Po 

Z y = ■ (28) 

Let us now build the variational state: 

|V>i,a) = K + ) - K">- (29) 
Eqs. (H) imply: 

% |Vi,a) = -e* ktru |Vi,a) 

Ti-Tr \lpl,A) = -Po IV'l, A) 

S* |^i, A ) = -o-l |^,a> (30) 

|V>i,a) has thus a wave- vector ki, a rotation quantum 
number p\ and a reflection quantum number a\ related 
to the quantum numbers of |?/>o) by the relations: 

ki = k + (tt, 0) 
Pi = -Po 

o-l = -o-l (31) 

It is thus a state orthogonal to the ground-state (even on 
a finite-size system where the topological sectors are not 
rigorously orthogonal). 

Since any local Hamiltonian has exponentially vanish- 
ing matrix elements between different sectors we have 

(tP+\H o \tPo)^0 (32) 

and | ipi t a > is thus degenerate with the absolute ground- 
state, their symmetries being related by relations (|3l|). 

We have computed numerically the spectrum of the 
QHCD model for three different triangular samples: N = 
20 = 4 x 5, N = 28 = 4 x 7 and N = 30 = 6 x 5 
sites. These samples are translation and TZ^ invariant 
(but without symmetry axis). As expected, the spectrum 
is exactly two-fold degenerate. Pairs of states satisfying 
Eq. ( |3l"| ) are exactly degenerate. The two ground-states, 
in particular, have k = (0, 0) 1Z n = 1 and k = (tt, 0) 

n n = -l. 
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3. Degeneracies in evenx even samples 

As remarked by Bonesteeli on a square lattice a n/2 
rotation exchanges sector (+, — ) and sector (— ,+) but 
sectors (— ,— ) and (+,+) remain inequivalent. A sim- 
ilar phenomenon occurs on the triangular lattice where 
27r/3 rotations permute cyclically 3 of the 4 sectors. Con- 
sider a finite-size triangular lattice with periodic bound- 
ary conditions which has the 27r/3-rotation symmetry 
(N = L x L for instance). One draws 3 cuts Ai which 
can be deduced from each other by ±27r/3 rotations. Of 
course only two of these three cuts are topologically dis- 
tinct but one can associate 3 parities to any short-rang 
dimer configurations anyway. These three cuts separate 
the sample in two sets of A^/2 sites. If, for instance, N/2 
is even the only allowed parities are (+, +, +), (+, — , — ) 
(—,+,—) and (—,—,+) and we recover 4 sectors as it 
should be. In this form it is now obvious that the three 
last sectors are symmetric by 27r/3 rotations and will 
have the same ground-state energy. One can also easily 
determine the relations between the quantum numbers of 
these triplets. As result, if N/2 is even, sectors (+,—), 
(— , +) and (— , — ) will be exactly degenerate and if iV/2 is 
odd (+, — ), (— , +) and (+, +) will be exactly degenerate 

If we consider how 27r/3 rotations and translations 
act on this three-dimensional manifold, one can eas- 
ily show that this 3-fold degeneracy can only take two 
forms, a) The degenerate states can correspond to a 
three-dimensional representation. The only possibility is 
that the corresponding momentum are Kj = i ( 2,3. b) The 
other possibility is that a two-dimensional representation 
is degenerate with a one-dimensional representation. The 
only possible realization is a k = R 
generate with two k = R 



2ir/3 



1 state de- 



'2tt/3 



j and j complex 



conjugate states. For these samples, the 27r/3 rotation 
plays exactly the same role as the one-step translation in 
the case of evenx odd samples. 

All these properties were checked numerically in the 
QHCD model on the triangular lattice. We found that 
the second case of degeneracy is realized for this model in 
the N = 28 sample (i.e. the four eigenstates of the quasi- 
degenerate ground-state multiplet have wave vector k = 
(0, 0) but belongs to different IRs of C3). In fact, it must 
be so for any evenx even sample at the RK point since the 
wave functions clearly have k = in each sector (equal 
amplitude superposition of all dimer coverings). This 
remains true as long as the gap does not close and should 
be verified in the whole spin-liquid phase of the model. 
In that case both Proposals B and B' of subsection II A 
are incorrect. 

We can ask what is the result of Oshikawa's insertion 
of a flux in such an evenx even spin-liquid. When bound- 
ary conditions are twisted, the Hamiltonian is no longer 
i?2ir/3- s y mm etric and the 2 irreducible representations in 
which the topological ground-states lie (-R27173 = 1 an d 
i? 2 Tr/3 = jij 2 ) merges into a single one with momentum 
k = (0,0). For this reason, during the adiabatic process, 
these four energy levels cannot cross and must necessarily 



map onto themselves after the complete flux quantum in- 
sertion. In that case Oshikawa's procedure do not relate 
the different ground-states with each other. 



B. 4-fold degeneracy in RVBSL phases 

The numerical data on MSE model suggest that the 
ground-state degeneracy is 4 in a spin liquid phaseElj, 
whatever may be the shape of the sample (Fig. ^). In 
the case of evenx odd sample [N — 30 = 6 x 5 in Ta- 
ble III), the quasi-degenerate levels exactly satisfy the 
constraint of Eq. (j3l|). In evenx even samples which have 
the 27r/3-rotation symmetry (N = 28 and 36 in Table III 
and Fig. 0), the quasi-degeneracy is observed between 
a k = (0,0) state and three k = Kj.,2,3, in perfect agree- 
ment with the arguments developed in paragraph IV A 3 . 
These are additional arguments in favor of the topolog- 
ical interpretation of the quasi-degeneracy observed in 
this modeled 

The same situation appears in the liquid phase of the 
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FIG. 10: First eigenstates of the multiple-spin exchange 
model on a 6x6 sites sample (parameters J2 = —2 and J4 = 1; 
the system is in a spin liquid phase). The first four levels with 
total spin S — and some S = 1 are displayed. The quan- 
tum numbers of the eigenstates are displayed at the right of 
the figure. The ground-state belongs to the trivial representa- 
tion of the group of translations, rotations and reflections (in 
particular k = (0, 0), po — 1 and <Tq = 1). The first S = ex- 
cited state is three- fold degenerate (wave- vectors Kj, p\ = —1 
and a\ = 1). The finite-size scaling indicates that this state 
collapses to the absolute ground-state in the thermodynamic 
limit (Table [o] and Fig. ^|). The third and fourth energies 
in the S — spin sector are also probably degenerate in the 
thermodynamic limit (global 4- fold degeneracy). 
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N 


28 


30 


32 


36 


k 


= (o, 0) 


-110.965 (1) 


-117.896 (-1) 


-126.608(1,1) 


-142.867 (1,1) 


ki 


= (tt,0) 


-109.396 (-1) 


-117.193 (1) 


-125.684(-1,-1) 


-142.615 (-1,1) 


k 


= (0,0) 


-108.168 (1) 


-116.590 (1) 




-142.098 (-1,1) 


ki 


= (7T,0) 


-107.633 (-1) 


-116.544 (-1) 




-141.883 (-1,-1) 



TABLE III: Energies and quantum numbers of the first two multiplets in the spin liquid spectra of the MSE model. The wave 
vector is indicated in the first column, the quantum numbers p* and a y are indicated in parenthesis after the energy. Samples 
N = 28 and 30 have no axial symmetry. Relations ( |3l| ) are satisfied in the even-odd . * stands for data where the 2 nd multiplet 
does not appear clearly, due to the absence of C3 symmetry of the sample. 



QHCD model on the triangular lattice (we consider the 
,7=1 and V — 0.95 point to avoid difficulties with the 
first order transition to the staggered phase at V = 1). 
This degeneracy property was altaady inferred by various 
authors at th e end of the eightiesBcrtl From the results of 
section [TV A we know that the ground-states inside each 
of the 4 topological sectors give rise to two energy lev- 
els (two doublets or one singlet plus a triplet depending 
on the sample symmetry). Indeed a quasi four- fold de- 
generacy of the ground-state is observed. More precisely, 
the two levels approach each other as the linear size is 
increased. For instance, this small level splitting drops 
from 5 = 0.02285 on a 16 = 4 x 4 sample to 0.008147 for 
30 = 6 x 5 sites. In samples without 2 it/ 3- rot at ion sym- 
metry a set of four nearly degenerate levels is observed 
(N = 24 = 6 x 4 for instance). 

We make the following assumption: a) the ground- 
state can be described in a short-ranged dimer basis, b) 
All n-dimcr correlations (n — 2, 3, • • •) are short-range 
and the corresponding correlation lengths are bounded, 
c) The Hamiltonian is local. From hypotheses a) and c) 
it is clear that for a large enough system the four topo- 
logical sectors are not mixed in the ground-state and the 
spectrum can be computed separately in each sector. Wc 
do not have any symmetry operation which connects all 
four sectors and we need a physical argument to explain 
the fact that energies should be the same in each sector 
(in the thermodynamic limit). Because of their topo- 
logical nature, it is not possible to determine to which 
sector a dimer configuration belongs by looking only at 
a finite area of the system. In other words, any dimer 
configuration defined over a large but finite part of the 
system can be equally realized in all sectors. The Hilbert 
space available to the system is the same over any fi- 
nite region of the system. In the absence of any form of 
long-range order the system can therefore optimize all its 
correlations with an arbitrary high accuracy equally well 
in each sector. At this point we can only conclude that 
the four sector will have the energy density and we can- 
not exclude the existence of a gap between the different 
topological sectors. However the numerical results ob- 
tained in the QHCD and MSE models indicate that it is 
not the case and that the four ground-state have asymp- 
totically the same total energy. We think that this should 
be true for a general short-range RVB spin liquid. Based 



on the relationship between the effect of a twist and the 
topological degeneracy in such systems, a vanishing gap 
between sectors in the thermodynamic limit is likely to 
be related to th e com plete absence of sensitivity to such 
a twist (section |III C ) . 



C. Miscellaneous remarks on the RVB 
ground-state degeneracy 

- Dimers and twist operator. The variational state 
IV-'i.a) can be deduced simply from | "0o) by changing the 
sign of the dimers crossing the cut A. Such an operation 
can also been seen as a 2tt twist of the spins of column 
0. As was noticed by BonesteelQ for the one-dimensional 
problem such an operation is mathematically related to 
action of the LSMA twist operator. From the physical 
point of view the reason why |^>i,a) has asymptotically 
the same energy as ji/'o) becomes clear: in the absence 
of stiffness and of long-range spin-spin correlations, the 
perturbation induced by the boundary term of Eq. ( |l4| ) 
cannot propagate and does not change the energy of the 
initial state: its only effect is to change the relative phases 
of the different topological components of the wave func- 
tion, and consequently the momentum and space symme- 
try quantum numbers of the initial state of the even-odd 
samples. 

- Fractionalization and topological degeneracy. To 
our knowledge all present theoretical descriptions of frac- 
tionalized— excitations in 2D magnets or related prob- 
lems^'tSESoEll) (we should also mention topological 
properties of Laugh linjs wave function for fractional 
quantum Hall effectO^j) imply topological ground-state 
degeneracies. In such pictures, the physical operation 
which transforms a ground-state into another is the vir- 
tual creation of a pair of spinons (by dimer breaking) 
followed by its annihilation after the circulation of one 
of them around the torus. In such a process a n phase- 
shift is introduced between the topological sectors (as in 
the above recipe). For samples with an odd number of 
rows this operation connects eigenstates with different 
k vectors (and space quantum numbers) as described in 
Eqs. (pl|). Oshikawa's adiabatic construction is of the 
same nature. 

- Numerical studies with a different topology. An in- 
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teresting check of the pure topological nature of this de- servable is: 
generacy could be obtained by studying the problem no 
more on a torus but on a surface with a different genus. 
On a sphere we expect an absence of degeneracy. Un- 
fortunately if a lattice can be represented on an infinite 
plane, both the number of links L and plaquettes P de- 
pend linearly on the number of sites N and Euler's rela- 
tion P — L + N = 2 — G constrains the genius G to be 
2 ! The torus is the only possible topology if we require a 
full translation imjariance in both directions. In a recent 
work, Ioffe et alEl have studied the absence of sensitiv- 
ity to disorder as an evidence for topological phenomenon 
in the liquid phase of the QHCD model on the triangu- 
lar lattice. They also used open boundary conditions to 
modify the topology of the system and argue in that case 
that the low-energy spectrum is free of edge states which 
could hide the actual ground-state degeneracy. 

— Example of RVB phase with 2 spins in the unit cell. 
A spin liquid state, seeming very similar to the state ob- 
served in the MSE model on the triangular lattice, has 
been observed in the J\ — J2 model on the hexagonal 
latticeE] for J\ = —1, J 2 = 0.3. No quasi-degeneracy 
of the ground-state has been noticed. It should be re- 
marked that in this system there are 2 spins \ per crys- 
tallographic unit cell and no degeneracy is expected on 
the basis of the topological arguments. 



D. Symmetry breaking in gapped phases 



N 

A = J2e iKl - T *Ps= (r j ,T j + u) (33) 

where Ps=o(rj,rj + u ) is the projector on the singlet 
state of two neighboring spins. A connects eigenstates 
with wave-vector ko to states with wave- vector ko + Ki 

On a finite size sample, with periodic boundary condi- 
tions, the expectation value of A is zero in any eigenstate, 
but < A 2 > could be non zero. If the order parameter V 
defined by: 



V 2 =<^g. S .\AU\^g. S . >/N 2 



(34) 



does not vanish in the thermodynamic limit, the system 
has columnar dimer long-range order with wave vector 
Ki. 



Let us now consider a perturbation of the Hamiltonian: 

H 5 = H -(8A + h.c). (35) 

At T=0, the intensive linear response on the observable 
A is measured by the susceptibility: 

X = If < \^ u 1 j* A^g-s. > (36) 



From the mathematical point of view, ground-state 
wave functions that break one-step translations or space 
group symmetries can be built from linear combinations 
of the degenerate ground-states of the even- odd sam- 
ples. In a completely equivalent way, ground-state wave 
functions that break rotation symmetry can be built in 
even x even samples. One could thus superficially con- 
clude that spontaneous symmetry breaking is possible in 
RVBSL, we will show below that this assumption is false. 

There are many features which show that this degener- 
acy property is a subtle one, both from the mathematical 
and physical viewpoints. 

The possible alternation of the spatial properties of 
the low-lying excitations with the parity of the number 
of rows of the sample (as observed in the QHCD model 
on the triangular lattice) is a first difficulty. The degen- 
eracy of the RVBSL is in fact quite different from that 
appearing in a VBC. We do not expect the VBC ground- 
state degeneracy to depend on the genus of the sample, 
as the RVBSL does. 

From the physical point of view also the two situations 
are quite different. An infinitesimal symmetry break- where 
ing perturbation is able to select one symmetry breaking 
ground-state of the VBC, but as we will show below this 
is impossible in the RVBSL. 

Let us call A the extensive non-diagonal observable 
appearing in the VBC in the thermodynamic limit. On 
a columnar VBC modulated in the u direction, this ob- 



This susceptibility is bounded from belo 
4P 4 N 2 



Jl: 



< X 



(37) 



where / is the oscillator strength: 



f = jj(1> B ...\[A,[H ,A]]\1> a . a .) (Msi 



The demonstration uses the properties of the spectral 
decomposition associated to the operator A: 

S(w) = -j- Y, l< ^a-s\An >| 2 S(u - w n ) (39) 



where uj n = E n - E g , s 



P 2 = i J S(co)du (40) 
Using the Cauchy Schwartz inequality one obtains: 

P 4 < ^ / uS(u)dLJ J u- l S{u)duj (41) 



ujS{uj)dui = f/2 



LO 1 S(ld)dld = x/2 



(42) 



(43) 



15 



which proves inequality (|37]). For a short-range Hamil- 
tonian the oscillator strength / is 0(1) and inequality 
( |37| ) implies that the T=0 susceptibility associated to 
a finite order parameter diverges at least as the square 
of the sample size: any infinitesimal symmetry breaking 
perturbation will select a symmetry breaking state. 

We will now show that for a RVBSL, where all the 
correlations functions are short-ranged with correlation 
lengths bounded by £, the susceptibilities of the medium 
remain finite in the thermodynamic limit. To do so we 
distinguish in Eq. |36| the contributions from the quasi- 
degenerate states of the topological multiplet (called 
| aj >) from the contribution of the other states of the 
spectrum, above the physical gap A. We thus obtain the 
following upper bound for the susceptibility: 



X = Xtop +XA 



(44) 



Xtop 



2 I I it \ a l >< a l\ At , , ,r-\ 

op =T7<^. S .l^ t T7 i ^A\i> g . s .> (45) 



N 



E — E 



XA < Jj£ [<i>g. S .\AU\lP g . s . > 

- <ip g . s .\A i \ai><a 1 \A\ip g . s .>] (46) 

where \cti > stands for the state(s) of the topological mul- 
tiplet connected to the absolute ground-state by A. Us- 
ing the local properties of A, A' \a\ > is in the same topo- 
logical sector as \a± > and < rjj g . s . \A' \a\ > is at most of 
0(N x 2- L / 2 ) (see paragraph [TV A lfr. As Eg. - E g . s . is 



supposed to decrease as exp{—L/£) (see Fig. |y|j x top goes 
to a constant when the size_Qf the sample goes to infinity, 
provided £ is small enoughEa. In a system with exponen- 
tially decreasing correlations, V 2 decreases as 1/N and 
Xa is trivially constant at the thermodynamic limit. 

In such a phase an infinitesimal field cannot induce a 
symmetry breaking and there could not be any sponta- 
neous symmetry breaking. 



V. CONCLUSION 

The main statements of this paper could be summa- 
rized in the following way. In any situations the ground- 
state of an SU(2) Hamiltonian with one spin-i per unit 
cell is degenerate in the thermodynamic limit. In fact 
the problem encompasses two very different cases: 

• Either the system has a T = macroscopic long- 
range order: in these situations the degeneracy of 
the ground-state encompasses and reveals the sym- 
metry breakings of the macroscopic ground-states. 
It is algebraic in the system size in case of a contin- 
uous SU (2) symmetry breaking and involves only 
a finite number of levels in case of a discrete sym- 
metry breaking. This situation is not restricted to 
spin-i Hamiltonians. 

• Or the system has no long-range order at T = 0, it 
is a spin liquid. In that case the degeneracy prop- 



erty depends essentially on the existence of an odd- 
half-integer total spin in the unit cell. In this paper 
we have only studied the case of a unique spin-i 
per unit cell, but in view of its topological origin 
we suspect the result to be true for the general case 
above. On the other hand, we have given one ex- 
ample with two spins \ in the unit cell and very 
plausibly a-upique ground-state in the thermody- 
namic limitEj. 

All these ideas had already been known by a number 
of authors. The really new contributions of this work are 
the following: 

• We have analyzed the Lieb Schultz Mattis Affleck 
conjecture for 2D spin systems and shown by the 
numerical study of a counter-example that a part 
of this conjecture is probably incorrect (section II). 

• For spin liquid systems, we have analyzed Oshikawa 
adiabatic construction and shown that it is equiv- 
alent to the hypothesis of an absence of stiffness of 
the gapped systems. From our numerical analysis, 
we conjecture that the spin liquids display in the 
thermodynamic limit a total absence of sensitivity 
to the boundary conditions (section |l|) . 

• We have shown that on a 2D torus with an odd 
number of rows, the ground-state is degenerate 
(which was already known), we have also demon- 
strated some exact relationship between the quan- 
tum numbers of these multiplets and extended the 
demonstration of degeneracy to some even x even 
systems with rotation symmetry (which is new). 
We have developed the idea that these spatial prop- 
erties of the spin liquid spectra depend crucially on 
the parity of the number of rows in the sample and 
have developed both exact and numerical examples 
to support this statement. The study of the QHCD 
model in even x even samples is an explicit example 
where proposals B and B' of LSMA and Oshikawa 
are not verified. 

• Endly we have shown why, in spite of the degener- 
acy of the ground-state, usual spin liquids cannot 
exhibit spontaneous symmetry breaking of any kind 
contrarily to VBCs (or chiral spin liquid) (the ulti- 
mate reason of our obscure and sometimes contro- 
versed previous claim that the invoked symmetry 
breakings in spin liquids were "nonphysical" ) . 

Acknowledgments: We acknowledge very fruitful 
discussions with S. Sachdev and V. Pasquier. One of us 
(C.L.) thanks the hospitality of I.T.P. (Santa Barbara) 
and the organizers of the Quantum Magnetism program. 
Special thanks are due to I. Affleck and D. Mattis whose 
vivid interest has prompted this work. This research was 
supported in part by the National Science Foundation 
under Grant No. PHY94-07194, and by the CNRS and 



16 



the Institut de Developpement des Recherches en Infor- 
matique Scientifique under contracts 010091-010076. Di- 
agonalization of the QHCD model were done on the Com- 
pacq alpha server of the CEA under project 550. 



when 



ki + k 2 = (ir [L y mod 2] 



CZ 

2tt^,0) 



If Sj? ot = this condition reduces to: 



(51) 



VI. APPENDIX A: SYMMETRIES OF 
SPECTRA WITH TWISTED BOUNDARY 
CONDITIONS 

A. Arbitrary twist angle (f> 

With periodic boundary conditions, spectra at momen- 
tum k and — k are identical because the Hamiltonian 
commutes with TZ^ the lattice rotation of angle ir about 
the origin. In the presence of a non zero twist, this rota- 
tion is no longer a symmetry : 



(47) 



If one defines the spin flip F by FS* = —S*F, then TZ V F 
commutes with and £^(k, S z ) — £^(-k, -S z ). 
For an arbitrary twist the only symmetries of are : 

• H„F 

• spin rotations about the z axis : e teStot 

• Translations 

• Other spatial symmetries which do not change the 
x coordinates. The lattice may have, for instance, 
an axis of symmetry parallel to the x direction. 

All these symmetry operators commute and the Hilbert 
space of Hrf, can thus be analyzed in terms of tensorial 
products of one-dimensional IRs. 



B. 



7T-tWlSt 



For 7r twist the spectrum has one additional symmetry. 
Since U(2mr)H^U(2mr)^ 1 — H^^nn for integer n , one 
finds: 

[U(2ir)U v ] H$ [U(2iT)n 7T r 1 = H 2 ^ (48) 
and therefore 



[U(2ic)K v ] H„ [[/(2tt)^] 



Hi, 



(49) 



This new symmetry induces extra degeneracies in the 
spectra. Using Eqs. (|4|) and (|l9|), one obtains: 



-E^=7r(Stot> k l) - E<t>=Tr(St ot , k 2 ) 



(50) 



ki 



(w, 0) L y odd 
(0, 0) L y even 



(52) 



Of course, these degeneracies are explicitly verified in all 
our numerical spectra. 

Finally, we address the issue of the property in the 
FTZtt transform of the eigenstates which are degenerate 
at <j) = (the quantum number associated to this prop- 
erty is rj = ±1). For this purpose we compute the com- 
mutation of FIZk with U(2ir) : 

F11 w U(2tt) = U(2Tr)FK„exp{2iiTSZ ot + 2iiTS£)(53) 
FTZ„U(27t) = U(2Tr)FK 7V (-l) L * Ly+Ly (54) 



where S% = J]p-n 1 #5- 



involves a sum over all the 



o — Z-*p=0 u n=a,p 

spins of the column 0. Since TZ^ trivially commutes with 
n we find that each doublet of degenerate eigenstates 
obey 7/i = {-l) L * L y +L yri2- 



VII. APPENDIX B: THE MULTIPLE-SPIN 
EXCHANGE MODEL AND THE TYPE I SPIN 
LIQUID PHASE 

The multiple spin exchange (MSE) Hamiltonian is an 
effective Hamiltonian which describes the exchange of lo- 
calized fermions on a lattice. In terms of multiple-spin 
permutation operators, the 2—, 3— and 4— spin exchanges 
can be reduced to : 



H 



j 4 J2 ( P4+p 4 



(55) 



where J| 



(resp. J4) measures a combination of the 
tunneling amplitude of 2— and 3— spin exchanges (resp. 
4— spin exchange). The permutation operators can be 
rewritten in terms of usual spin operators: Py = 2S^ • 
Sj + 1/2. P1234 + h.c. is a polynomial of degree two in 
Si • Sj . A complete expression in terms of spin operators 
can be found in Ref.Lj. This Hamiltonian has a spin- 
gap phase in a wide range of ferromagnetic J| , and 
antiferromagnetic J4 couplingscfl. The spin liquid spec- 
tra studied in the present work are at the specific point 



X 



eff 



-2 and J4 = 1, where the spin gap is of order 1. 
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The existence of a VBC phase for such a coupling is rather 
well docume nted 1 1 bi fh the exact nature of the order is still 
under debatdaMfeafei 

This softening due to quantum fluctuations always notice- 
able on the stiffness and spin-wave velocity, seems even 
larger on the maxima of the dispersion curves, as it has 
also be noticed on the square lattice£3. 
The first excited state of the spectrum is an S — 1 
k = (n, 7r) eigenstate: its gap to the ground-state goes 
to zero linearly with the size of the sample; its value is 
0.630 (resp. 0.575) in the N=32 (resp. N=36) sample. The 
first eigenstate with a wave- vector k = (n, 0) has a gap to 
the ground-state of 4.885 (resp. 4.850) in these two sam- 
ples: it is triplet (S — 1). The first eigenstate with S = 
and k = (ir, 0) has a huge gap to the ground-state of 5.195 
(resp. 5.075). 

This relation has also been used to generalize to S'tot 
the LSMA argument iiv-arder to discuss magnetization 
plateaus in one dimensionEfl. 

In fact we will show below that the four-fold degeneracy 
of the RVBSL can sometimes be realized in even x odd and 
even x even samples with different irreducible representa- 
tions of the space group. 

In small frustrated samples the absolute minimum could be 
observed for a very small, non zero, twist. This peculiarity 
is no more present for the stablest symmetric samples of 
size larger or equal to 28. 

Another equivalent definition is based on winding mimbers 
of transition graphs with a reference configuration!! 
A ground-ptiate degeneracy has been evoked by Wen and 
co-workersH in chiral spin liquids. To our understanding, 
this is a bit more complex than the present case. A chiral 
spin liquid is a gapped phase with long-range order in chi- 
rality (local observable defined as the triple product of 3 
neighboring spins). It not only breaks reflection symmetry 
but also time reversal symmetry and it is supposed to be 
2{k) d degenerate, where k is an even index related to the 
fractional statistics of the elementary excitations. 
Strictl y sp eaking £ should be < log 2 but going from Eq. ^l] 
to Eq. p2j dramatically overestimates the scalar product in 
the case of an RVBSL. The reason is that if two dimer 
coverings c + and c~ maximize ^c + |c~^ they only differ 
along a single large (~ L) loop. They have different local 
correlations along the loop and their energy difference is 
of order L and it is very unlikely that their weights in the 
ground states if> + (c + ) and ip~(c~) are both of order one. 



